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ABSTRACT 

Context. The relation between Galactic cosmic-ray electrons, magnetic fields and synchrotron radiation. 

Aims. We exploit synchrotron radiation to constrain the low-energy interstellar electron spectrum, using various radio surveys and 
connecting with electron data from Fermi-LAT and other experiments. 

Methods. The GALPROP programme for cosmic-ray propagation, gamma-ray and synchrotron radiation is used. Secondary electrons 
and positrons are included. Propagation models based on cosmic-ray and gamma-ray data are tested against synchrotron data from 22 
MHz to 94 GHz. 

Results. The synchrotron data confirm the need for a low-energy break in the cosmic-ray electron injection spectrum. The interstellar 
spectrum below a few GeV has to be lower than standard models predict, and this suggests less solar modulation than usually assumed. 
Reacceleration models are more difficult to reconcile with the synchrotron constraints. We show that secondary leptons are important 
for the interpretation of synchrotron emission. We also consider a cosmic-ray propagation origin for the low-energy break. 
Conclusions. Exploiting the complementary information on cosmic rays and synchrotron gives unique and essential constraints on 
electrons, and has implications for gamma rays. This connection is especially relevant now in view of the ongoing PLANCK and 
Fermi missions. 
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by the PAMELA experiment (lAdriani et alJ l2009h is also dis- 
cussed in these papers. Our preliminary study of the synchrotr on 
spectrum using GALPROP was given in lOrlando et al.l d2009t) . 

The main objective of the present paper is to constrain the 
interstellar electron spectrum using a combination of the lat- 
est electron spectrum measurements and synchrotron radiation, 
leaving the question of its origin via injection and propagation 
as a side-issue. Hence while we use the CR propagation code 
GALPROP to generate interstellar spectra for various propaga- 
tion scenarios, the latter are not our main focus. Nevertheless it 
is shown that some current models are actually excluded on the 
basis of the synchrotron data. The relation between synchrotron 
and electrons is complicated by the presence of secondary elec- 
trons and positrons, and hence these are included in our model 
and their effect is addressed. 

In previous analyses we were able to constrain the total B- 
field Bio, on the basis of synchrotron data and the distribution 
of cosm ic-ray electrons derived from gamma rays (IStrong et al.1 
2004a); we obtained a local value B t0 , - 6/iG and a scale length 
of 8 kpc in R and 1 kpc in z. The value of B tot has frequently been 
quoted in the literature as an independent measurement of the to- 
tal B-field. Since then the CR source distribution has been mod- 
ified to better reflect the distribution of SNR as traced by pul- 
sars (iStrong et alJl2004bl) . and this influences the derived radial 
variation of B tol ; the larger radial variation of the source func- 
tion implies a smaller variation of B tol . Also the new measure- 



1. Introduction 

Direct measurements of cosmic-ray (CR) electrons extend from 
TeV down to 1 GeV (and lower energies from spacecraft 
like Ulysses and Voyager), but solar modulation complicates 
their intepretation at energies below about 10 GeV. Fermi-LAT 
( Abd o et alj|2009bt lAckermann et al.l 12010) has made the cur- 
rently most precise electron measurements in the range 7 GeV - 
1 TeV; here modulation is of less importance although still sig- 
nificant in the lower part of this range. Synchrotron radiation 
(from tens of MHz to tens of GHz) probes interstellar electrons 
from 0.5 to 20 GeV for the typical Galactic magnetic field (here- 
after B-field) of a few fiG, and hence can be used in conjunction 
with direct measurements to construct the full spectrum from 
GeV to TeV. At low energies this will finally allow an indepen- 
dent estimate of solar modulation for testing heliopheric propa- 
gation models. In contrast such a probe of the low-energy inter- 
stellar spectrum is not available for CR nuclei. 

An extensiv e review of CR pr opagation including electrons 
can be found in Strong et al. (20 07|); a recent glo bal viewpoint 
for the Milky Way is giv en in | Strong et all (1201 Oh . and a propa- 
gation parameter study in lTrotta et alJ d201 ll) . Further discussion 
o f the electron spectrum measured by Fermi - LAT can be found 
inlDelahave et alJd2010l) ; lGrasso et al](l201 lh ; ldTBernardo et alJ 



(201 1). Its relation to the high-energy positron excess discoved 
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ment of the electron spe ctrum by Fermi -LAT (lAbdo et al . 2009b; 
lAcker mann et al. 2010) and the elimi nation by Fermi-L AT of 
the EGRET gamma-ray 'GeV excess' (lAbdo et aU l^Ogafl lead 
to an upward revision of the value of B lol . We concentrate on the 
spectral aspects of synchrotron emission in this paper. 

The magnitude of the B-field is a free parameter in this anal- 
ysis, since while the regular component can be determined from 
rotation measures of pulsars and extragalactic sources, this is 
only a fraction of the total field. Our approach is to use the mod- 
els for regular component derived from RMs, combining these 
with a random field to be determined, and the latter is one of the 
results of our analysis. 

Extensions to lower frequencies whe re absorption is im- 
portan t ha ve been made, see for exa mple IStrong & Wolfendalel 
dl97l and lWebber & Higbiel J2008h . 

In a paper complementary to this one, we also address the 
relation between cosmic -ray electrons, synchrotron and B-fields 

3 



2. GALPROP development 

A description of the GALPROF0 model can be found 
in IStrong et al.l (2007b an d refer e nces therein; in part icular 
see IStrong & M oskalenko ( 1998b, IStrong et al.l (|2004a) and 
iPorter et al Jd2008l) . and the GALPROP Explanatory Supplement 
available from the GALPROP website. 

A shortcoming of GALPROP up to now has been the sim- 
plified B-field modelling, using only a random component and 
a simple 2D exponential dependence. With the availability of 
excellent radio continuum surveys from tens of MHz to tens of 
GHz, including the WMAP satellite data, more sophistication is 
desirable; we have therefore introduced full 3D models for both 
regular and random B-field . For the present spectral study only 
the total field is required (derived from the 3D model), and a 2D 
propagation scheme is sufficient. 



2.1.1. Regular field 

The synchrotron emissivity (in erg s _! Hz 1 ) of an isotropic 
distribution of monoenergetic relativistic particles in a uniform 
magnetic field has polarized components parallel and perpen- 
dicular t o the projection of the field on the line-of-sight to the 
observer (Longaii-2010): 



V3 e 1 

e P ar(y) = -= jB perp [F(x) - G(x)] 

V3 e 3 

£perp(v) = — zBperp [F(x) + G(x)] 



(1) 



(2) 



where x = v/v c , with v c = -^B perp y 2 and with y the electron 
Lorentz factor, B perp = B(x, y, z) sina, with a the angle between 
the magnetic field and the line-of-sight. The fun ctions F(x) an d 
G(x) are defined in terms of Bessel functions (Longair 2010), 
with: 



Xinf 
Ks/i(x') dx' 

G(x) = x K 2 /i(x) 



(3) 



where Ks/3(x) and Ky^ix) are the modified Bessel functions 
of order 5/3 and 2/3. They are conveniently provided as C li- 
brary functions in the GNU Scientific LibrarjQ. The resulting 
synchrotron spectrum has a broad maximum centred roughly at 
the frequency v c and the maximum has a value v mflx =0.29 v c 
(lLongairl2010h . 

The polarization formulation will be of use for our future 
work. Here we are only interested in the total intensity given by 
the sum of the two components described above: 



s(v,y)= V3 -^B perp F(x) 



(4) 



2.1. Synchrotron emissivity calculations 

From the spectrum of particles (here electrons or positrons) com- 
puted by GALPROP at all points on the 3D grid, we integrate 
over particle energy to get the synchrotron emissivity for the reg- 
ular and random fields. 

The emissivity as seen by an observer at the solar position is 
computed as a function of (x, y, z, v). The spectrum and distri- 
bution of the emissivity thus depends on the form of the regular 
and random components of the magnetic field, and the spectrum 
and distribution of CR leptons. 



1 our earlier work was based on an electron spectrum adjusted to fit 
this excess, and which was a factor 4 higher than locally measured, 
while now we use the measured electron spectrum. 

2 iJaffe et all i201lh . Both papers use the GALPROP code to model 
CR propagation. The main difference is that they consider emission in 
the Galactic plane from 408 MHz to 23 GHz, while we use data out 
of the plane down to 22 MHz. They address the spatial variations in- 
duced by the B-field, while we concentrate only on the spectral aspects. 
We explicitly compare the contributions from secondary leptons, while 
they consider only the total. They use polarization (with Faraday rota- 
tion) to separate the regular and random B-field components. They also 
use rotation measures to constrain the regular B-field. Despite the dif- 
ferences in approach, the conclusions - on issues common to the papers 
- are consistent. A rec ent paper describing a similar approach to ours 
(Bringmann et al. 201 1) has been brought to our attention. 

3 The code is publicly available at http://galprop.stanford.edu 



2.1.2. Random field 

For a randomly oriented field the emissivity is isotropic and ob- 
tained by integrating the regu lar field expressions o ver all solid 
angles. The result is given by dGhisellini et al.l fl988) 

£rand( v ) = C X 2 [K4/iK l /j, - -*(-^4/3&*/3 - K1/3K1/3)] (5) 

x = v/Ve, v c = ^-t c B ran y 2 , C = 2V3^rB ra „ erg s" 1 Hz" 1 , 
and the Bessel functions £4/3, ^1/3 are again computed using 
the GNU Scientific Library. 

Our implementation has been checked by integrating the reg- 
ular field expression over solid angle, giving exact agreement 
with this formula. 

2.2. Synchrotron intensity 

With GALPROP calculation of emissivity on the grid, we inte- 
grate over the line-of-sight to get the synchrotron intensity for 
the regular and random fields. The synchrotron intensity at fre- 
quency v is then given by 

/(v) = f e(v) ds (6) 



4 http://www.gnu.org/software/gsl 
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The observed brightness temperatur^ of the radiation seen in a 
given direction is 



T(v) oc 



c 2 /(v) 
2 v 2 



(7) 



The resulting synchrotron skymaps for a user-defined grid 
of frequencies are output by GALPROP in Galacti c coordinates 
either as CAR (Carre projection) or in HEALPix (Gorsk iet al.1 
120051) (the preferred format). The emissivity as seen by an ob- 
server at the solar position is also output as a function of (R, z, v) 
(spatial 2D) or (x,y,z, v) (spatial 3D). 

2.3. Galactic magnetic field models 

In the same way as for the other Galactic constituents (gas, ISRF, 
cosmic rays), the magnetic field is defined on a grid, in 2D or 3D. 
Only the 3D case is relevant for the full B-field model, although 
the 2 D case is retained for compatibility. 

In Strong et al.l d2000h only the random component of the 
magnetic field was present and was implemented in 2D with an 
exponential law for the component. 

Since that work, many 3D models of the Galactic magnetic 
field have been implemented in GALPROP in order to calculate 
the synchrotron emission from the Galaxy. The regu lar B-field 
used i n the present work is the model RING-ASS of ISun et"al] 
(2008) for the disk, based on rotation measures of extragalactic 
radio sources. This has typically B re „ — 2pG. A toroida l halo 
field is also included as prescribed in Sun & Reich! (1201 Of) , hav- 
ing a typical value of 2fiG. We include the regular field in order 
to make our model compatible with current information, but in 
fact since it is much less than the random field, this is not critical 
to our study. 

From Section 2,v c = (E/l GeV) 2 (5 ra „/7.5yuG)x 240 MHz, 
so that the full range of synchrotron frequencies used, 22 MHz 
to 94 MHz, traces electron energies from roughly 0.5 to 20 GeV 
for our adopted B-field. 

At 408 MHz and above, secondary leptons (above a few 
GeV) become less important for synchrotron and the relevant 
leptons are measured directly without much solar modulation, 
while at higher frequencies (>1 GHz) free-free emission can 
start to enter, so this is the best frequency to determine the ran- 
dom B-field. Our model for B nm is therefore based on fits to the 
data at 408 MHz. The random field is modelled as a double expo- 
nential in (R,z), the free parameters being the two scale lengths 
(30 kpc in R and 4 kpc in z) and the local B-field: B mn = 7.5yuG@. 
This model reproduces the longitude and latitude distribution of 
synchroton at 408 MHz sufficiently well for our purpose. Note 
that since we are mainly concerned with spectral shape in this 
paper, the absolute value of the B-field and its spatial distribu- 
tion are not critical, and affect mainly the relation of electron 
energy to synchrotron frequency. 



3. Cosmic-ray model 

We use the GALPROP models described in lStrong et ail d2010h . 
to which we refer for details. These models have been adjusted 
so that the propagation parameters are consistent with CR nuclei 



5 In this paper we use the temperature spectral index, as is commonly 
used in radioastronomy (e.g. literature in Appendix A), which is equal 
to to the intensity spectral index plus 2. 

6 We note that lSun et "all d2008h use an electron spectrum a factor 3 
higher than the Fermi-LAT measurements, and hence obtain a random 
field lower than ours (3/iG). 



secondary-to-primary ratios. Only the electron injection spec- 
trum and the B-field are varied with respect to these models; 
these do not affect the validity of the propagation parameters. 

Electrons and positrons lose energy by synchrotron radia- 
tion, and this is included in GALPROP self-consistently us- 
ing the B, ot of the adopted model. Energy loss by inverse 
Compton scattering, ionization, Coulomb and bremsstrahlung 
are also included, although the latter three processes are of mi- 
nor importance at the electr on energies of interest here (see 
[Strong & Moskalenkc] d!998l) for a plot of the loss processes). 

We concentrate on the spectrum of electrons as determined 
by synchrotron data, combining these with direct measurements. 
Electron spectra are modelled to be consistent with Fermi-LAT 
data in the range 7 GeV to 1 TeV, including an estimate of 
solar modulation which still has some effect at these energies. 
Below 7 GeV the primary electron injection spectrum is mod- 
elled with 2 breaks. The minimal constraint is that the propa- 
gated spectrum must not be below the directly-measured spec- 
trum (i.e. at least some solar modulation is present). Secondary 
electrons and positrons (from pp, p-He, He-p and He-He interac- 
tions via pion-decay) are included using the standard GALPROP 
treatment with locally-measured proton and Helium CR spectra 
as a normalization. Note that since Fermi-LAT measures elec- 
trons plus positrons, our primary electrons source effectively in- 
cludes the hig h-energy positrons p rimary component measured 
by PAMELA ( Adriani et al. 2009), so the lat t er is not explicitly 
model led here (see e.g. iGrasso et a i Bernardo et al. 

(1201 lh for such models). 

As usual all the GALPROP parameter files for the models 
will be available as Supplementary Material. The plots will be 
made available in numerical form on request. 



4. Radio data 

4.1. Review of information on spectral indices 

A revie w of radio con t inuum surveys and their calibration is 
given in lReich & Reichl d2009l) .We are interested in both spectral 
and spatial properties of the synchrotron sky, but not in the fine 
angular details. For the spatia l distribution the most useful is the 
full-sky 408 MHz survey (lHaslam et al.lll982h . since it has full 
sky coverage, a well-established calibration and zero level, and is 
a standard in this context. It has also the advantage that contribu- 
tions from non-synchrotron components (eg free-free emission) 
are rather small. For spectral information we have assembled a 
set of surveys which are described in Section 4.2. 

There are also many detailed observational studies aimed at 
measuring accurate sky temperatures and spectral indices at par- 
ticular wavelengths, for particular sky regions. These are very 
valuable in providing absolute values, especially at high Galactic 
latitudes. They also address the question of the contribution from 
extragalactic radio sources. We list some representative results 
from the literature on spectral indices in Appendix A. In sum- 
mary, there is a wealth of measurements showing that the spec- 
tral index of the synchrotron emission increases steadily from 
about 2.5 to 3.0 over the frequency range from tens of MHz to 
tens of GHz. While there is considerable scatter in the actual 
values, they still provide an essential observational constraint on 
any model for the CR electron + positron spectrum. Here we 
use a representative sample to compare with our models, but the 
comparison is just indicative of the general trend since the ex- 
periments cover many different sky areas. 
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4.2. Radio surveys used 

We use surveys at frequencies from 22 MHz to 23 GHz (and up 
to 94 GHz for WMAP) to compare directly the synchrotron spec- 
trum with the models. The surveys used are summarized below. 
Some were obtained directly from their authors, others from the 
Bonrfl and LAMBDA0 websites. The combined zero level and 
extragalactic/CMB corrections were taken from the literature as 
stated below. 

22 MHz: DRAO Northern hemisphere s ur vey: Roger et al.l 
(fl999h. 45 MHz: North: M aedaetalJ (119991). South: 



Alva rez et al.l (Il997l) . combined all skv: iGuzman et al.l (120111) 
This is complete apart from two regions of 10° and 20° 
radius out of the plane. An offset of 550K was subtracted 
dGuzmanetal.1 120111). 150 MHz: Parkes- Jodrell Bank all sky 
survey: lLandecker & Wielebinskil d!970l). 408 MHz : Bonn - 
Jodrell Bank-Parkes all sky surve y: [Haslam et alj (fl 982). 
An offset of 3.7K was subtracted dReich & Reichl 1 198 8 a PI. 
1420 MHz: Stackeri-V i lla El iza all s ky survey. NorthT[R eich 
dl982l) : iReich & Reich] d 19861) South: iReich et alj 4200JJ). An 
offset of 2.8K was subtracted dReich et alJl2004 2326 MHz: 
Rhodes southern hemisphere survey: Jonas et al] (1 1998b . This 
was not used here due to restricted coverage, but will be used 
when the analysis is extended to lower Galactic latitudes in 
future. 23 - 94 GHz: For WMAP we used the spectral-index 
maps generated on the basis of WMAP polarized data by 
Miville-Desch enes et al.l (120081) , which we used to scale the 
408 MHz map to the WMAP frequencies^ in this way the 
synchrotron radiation is extracted essentially uncontaminated 
by thermal and spinning dusQ- 

The surveys used have full-sky (or almost full-sky) coverage 
apart from 22 MHz, but here the coverage is still almost com- 
plete in the sky regio n used. The coverage of these surveys can 
be seen graphically in lde Oliveira-Costa et al.l d 20 08 r^l 

In order to avoid absorption effects at low frequency, and 
free-free emission at higher frequencies, and to avoid effects 
of zero-level corrections and local emission, this analysis is 
restricted to regions out of the Galactic plane but avoiding 



7 http://www.mpifr-bonn.mpg.de/survey.html 

8 http://lambda.gsfc.nasa.gov 

9 IReich et al.l d2004l) give 2.7K. lGuzman et al.1 d2011l) give 1.6K, but 
the difference is not significant here. 

10 iMiville-Deschenes et alj d2008h used the 408 MHz maps uncor- 
rected for zero offset (Miville-Deschenes, private communication), so 
we use the same maps for consistency when applying their index maps. 

11 Their re-analysis of the WMAP synchrotron data including spin- 
ning dust correction resulted in a lower intensity than previous analysis. 
Spinning dust emission is produced by small grains rotating and pro- 
duces unpolarized radio emission. They analyzed the combination of 
the WMAP polarization and intensity data finding strong evidence for 
the presence of unpolarized spinning dust emission in the 20-60 GHz 
range. They performed an analysis of the WMAP synchrotron emis- 
sion at 23 GHz where the signal to noise ratio is the highest and the 
polarised emission is only synchrotron. Their estimates of the intensity 
at this frequency are based on extrapolation of the Haslam 408 MHz 
data with a spatially varying spectral index constrained by the WMAP 
23 GHz polarization data. Hence, supposing that the synchrotron spec- 
tral index does not vary with frequency over the WMAP range, they 
found an anomalous emission with a spectrum from 23 to 61 GHz in 
accordance with the models of spinning dust. 

12 They give an extensive list of radio surveys including older ones. 
They also provide some of them in HEALPix and use a subset for a 
principal components analysis with the aim of predicting the radio sky 
at any frequency. However issues of calibration and zero-level are not 
addressed by them so we do not use their results here. 



the polar regions, specifically 10° < \b\ < 45 °. The thermal 
contribution even at 1420 MH z is then small dReich & Reichl 
1988b; Broad bent et all 1 19891) . a recent estimate being 15% 
dDickinson et al.ll2003h . and less at lower frequencies (while at 
WMAP frequencies our maps separate out the non-thermal com- 
ponent as explained above). We also avoid the North Polar Spur 
by avoiding the region 340° < I < 40°, although this hardly af- 
fects the result over such large sky regions, as we have verified. 

It is worth noting that lower frequency dat a are available (e.g. 
down to 1.3 MHz f rom the RAE2 satellite. iNovaco & Brownl 
d 19781) ). discussed in lStrong & Wolfendalel d 19781) . but these are 
strongly affected by absorption so are not used here. 

For the spectral comparisons, the s urveys and model s pre- 
dictions were converted to HEALPix dGorski et al.l I2005T) and 
averaged over the stated sky region. Since HEALPix has equal 
solid angle pixels the correct averaging is ensured. 



5. Electron and positron data 

The CR electron and posit r on da ta used are as foll ows: AMSOl 
dAMS Collaboration et aD 120021). CA PRICE94 dBoezio et al 



2000), HEAT (DuVernois et al. 2001), SANRIKU (Kobavashi 
M999h . BETS dTorii et al.1 120011). PPT-BET S dYoshida et al 
2008). ATIC-1-2 dChang et al.ll2008l). H.E.S.S. dAharonian et al 



2008, 120091) . Fermi-LAT (Abdo et al. 2009b; Ackerman n et al 



2010), PAMELA dAdriani et al.ll201 ll). The data are taken from 
the CR database described in [Stron g & Mos kalenkold2009l) R 

Fermi-LAT is taken as definitive above 7 GeV, while the ex- 
tension to lower energies is based on AMSOl, CAPRICE94 and 
HEAT. The low-energy data were taken around solar minimum 
with corresponding modulation levels. For HEAT the modula- 
tion levels quoted were 755 MV for 1994 and 670 MV for 1995 
dDuVernois et al.l|2001l). The C APRICE flight was in 1994, and 
quoted 600 MV dBoezio et al.ll2000t) . The AMSOl flight was in 
1998 at the end of a solar minim um, and quoted 650 + 40 MV 
dAMS Collaboration et ai1l2002l) Q The Fermi-LAT data was 
taken in 2008-2010 during an extreme solar minimum of polarity 
opposite to that of 1994-5, but for which no reliable modulation 
level is available. The PAMELA data were taken during 2006- 
2010, and a modulation of 600 MV is used in their paper. 

Fermi-LAT measured the electron spectrum from 7 GeV to 
1 TeV, with unprecedented accuracy. Above 20 GeV the Fermi- 
LAT electron spectrum can be fitted with a simple power law 
with spectral index 3.04. H.E.S.S. measured the electron spec- 
trum from 340 GeV to 5 TeV with a tendency to be higher than 
Fermi/LAT measurements in the region of overlap but in agree- 
ment within the systematics, and with a steepening above 1 TeV. 
ATIC - 1 - 2 measured the spectrum from 20 GeV to 2 TeV and is 
consistent with Fermi measurements up to 300 GeV, but found 
a peak between 300 and 700 GeV, which is not confirmed by 
Fermi/LAT. Energies above 100 GeV are not important for syn- 
chrotron (the corresponding frequencies lie beyond 1 THz ) but 
are included in the model and data for completeness. 

The el ectron spectrum (1- 625 GeV) recently measured by 
PAMELA dAdriani et al.ll2.Ql ll) is generally consistent with that 
measured by Fermi-LAT, although it is about 20% higher in the 
7-20 GeV range of energy overlap. The overall spectrum is also 
slightly softer than found by Fermi-LAT (and consistent with a 
rising positron fraction, see discussion in Adria nTet al.1 d201lh ). 
We choose here to baseline on Fermi-LAT, while bearing in 



13 available from http://www.mpe.mpg.de/~aws/propagate.html 

14 The modulation levels quoted by the authors implicitly assume an 
interstellar CR spectrum, so the derivation may be a circular argument. 
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mind this difference, which does not affect our conclusions. We 
do not address t he high-energ y positron excess reported by the 
PAMELA group lAdriani et aLl d2009h . since the main excess lies 
above the energies of importance for synchrotron, and in addi- 
tion because the absolute PAMELA positron spectrum is not yet 
available. 

Lower energy (100 MeV and below) electron data are avail- 
able from the Voy ager (I Webber & Hi gbie 2008) and Ulysses 
dHeber et al.l 12005) spacecraft, but they are not constrained by 
synchrotron and will not be used here. 



6. Results 

6. 1 . Pure Diffusion Model 

We consider first a 'pure diffusion' model, with a halo height of 
4 kpc. The comple te set of GALPROP parameters are given in 
IStrong et ail (1201 Ol) . model z04LMPDS. The electron injection 
spectrum breaks at 4 GeV and 50 GeV, with indices 1.6/2.5/2.2. 
The break at 50 GeV is to fit the Fermi-LAT low-energy upturn, 
the break at 4 GeV to fit low-frequency synchrotron. A cutoff 
at 2 TeV is introduced to reproduce the H.E.S.S. data, although 
this has no effect for the synchrotron, the corresponding frequen- 
cies being far too high. Fig Q] shows the interstellar electron and 
positron spectra for this model, and also for various modulation 
levels using the force-field approximation. Fig [2] shows that this 
model gives a reasonable fit to the synchrotron spectrum, and il- 
lustrates the role played by the secondary leptons which have a 
steeper spectrum than primaries and contribute significantly to 
the low-frequency synchrotron. Fig [3] shows the spectral index 
as a function of frequency for this model, for primary and sec- 
ondary leptons and for total leptons. 

We next show the effect of varying the low-energy (<4 GeV) 
electron injection index, from 1.0 to 2.5. Fig [4] shows the in- 
terstellar electron spectra for these models, and also for various 
modulation levels using the force-field approximation. It is clear 
that the electron data alone cannot distinguish the models due to 
the modulation uncertainty, so that the synchrotron constraints 
are essential. 

Fig|5]shows that a low-energy primary electron injection in- 
dex of 2.0 is at the limit of the low-frequency synchrotron data. 
The best fit is actually for an injection index around 1.3. Fig [6] 
shows the synchrotron spectral index for these models, com- 
pared to values from the literature described in Section 4.1. 

In Fig [7] the primary electron spectrum has been cut off be- 
low 4 GeV to illustrate the contribution from those energies; 
since removing these low-energy electrons eliminates most of 
the low-frequency synchrotron from primaries, it shows that 
low-frequencies (below 100 MHz) are dominated by leptons 
with energies less than 4 GeV. Secondary leptons produce one 
third of the observed low-frequency intensity and hence make 
determination of the primary spectrum more difficult. Secondary 
leptons together with primaries above 4 GeV already account for 
50% of the low-frequency synchrotron. 

6.2. Reacceleration model 

We now consider a reaccleration model, also with halo height 
4 kpc. The comple te set of GALPROP parameters are given in 
IStrong et al.l d2010l) . model z04LMS; the injection spectral in- 
dex above 4 GeV has been reduced from 2.42 in that model 
to 2.3 to better fit the Fermi electron data above 20 GeV. The 
range 7-20 GeV shows a slight steepening in Fermi-LAT and 
PAMELA data, but no attempt has been made to reproduce this 
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Fig. 1. Electron (upper) and positron (lower) spectra for pure 
diffusion model with primary low-energy electron injection in- 
dex 1.6. Modulation 0=0,200,400,600,800 MV. NB Fermi-LAT 
includes positrons. Cyan open circles: AMS01; green crosses 
and filled circles: CAPRICE; blue squares: HEAT; red filled cir- 
cles: Fermi-LAT; black filled circles: PAMELA; blue triangles: 
SANRIKU; red crosses: BETS, PPT-BETS; cyan open circles: 
ATIC-1-2; green filled and open squares: H.E.S.S. For references 
see text. 

in the reaccelation model since we want to test an existing pub- 
lished model; it has no effect on our conclusions. As in the case 
of the pure diffusion model a cutoff above 2 TeV has been in- 
troduced to fit the H.E.S.S. data, although this has no signifi- 
cance for synchrotron. The lepton and synchrotron spectra for 
this model are shown in Fig[8] Fig[9]and the synchrotron indices 
inFigJTU] 

It is clear that this particular reacceleration model is not con- 
sistent with the observed synchrotron spectrum, since the sum 
of primary and secondary leptons produces too high intensi- 
ties at low frequencies, and the low-frequency spectral index is 
too large. It could be adjusted by making the low-energy injec- 
tion index smaller, as for the pure diffusion model. However a 
large part of the excess comes from the secondary leptons which 
have a large peak due to reacceleration which makes them equal 
to primary electrons around 1 GeV, and which cannot be ad- 
justed very much in this model; this peak is not present in the 
pure diffusion model (see comparison for secondary leptons in 
Fig 0. Decreasing the B-field can improve the low-frequency 
fit but then the high-frequencies are under-predicted since the 
overall spectrum is steeper than the model predicts. Only if the 
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Fig. 2. Synchrotron spectra for pure diffusion model with pri- 
mary low-energy electron injection index 1.6. Synchrotron from 
primary electrons (upper), secondary leptons (middle) and total 
(lower). For synchrotron data references see Section 4.2. 



secondaries are removed does the synchrotron give a good fit, 
while the secondary production is certainly present. Arguing dif- 
ferently, we note that the secondaries already produce the low- 
energy synchrotron intensities, which would preclude the exis- 
tence of electron primaries. Either way the model is problematic 
for synchrotron. This does not mean that reacceleration models 
are excluded by this study, but it does pose a challenge for fu- 
ture work on such models. Reacceleration is surely present at 



some level on physical grounds, but probably less than currently 
adopted on the basis of B/C data. 

This is an example where solar modulation can be invoked 
to get agreement of a model with directly-measured CR data (in- 
creasing the interstellar spectrum and modulation appropriately), 
but not for synchrotron which probes the interstellar spectrum. 

6.3. Low-energy spectral index from propagation ? 



The low-energy injection indices deduced for the standard prop- 
agation models are unexpected (see Discussion), so it is valid to 
ask whether they could be produced by propagation using a more 
'normal' injection spectrum with index 2. One way to do this is 
to reduce the energy losses by making the propagation region 
smaller - which means reducing the halo size. Then the steepen- 
ing by propagation is reduced. This is at the expense of no more 
fitting B/C, 10 Be/ 9 Be and local electron measurements, but it is 
worthwhile illustrating this explicitly. FigQT]shows a model with 
halo height 1 kpc, but otherwise the same as the previous pure 
diffusion model (with its low-energy injection index 2 and halo 
height 4 kpc). (The diffusion coefficient has not been changed to 
fit B/C since this would increase the propagation time and the 
energy losses would be unchanged from the larger halo case, 
and the synchrotron spectrum would be the same as before.) 
To get the correct synchrotron intensity, the electron flux has 
to be increased by 3 compared to that observed locally to com- 
pensate the smaller integration length (increasing B would just 
lead back to the same losses again). With this unnatural scenario 
we can indeed reproduce the observed synchrotron spectrum 
(FigfTTTi. However since several other constraints are thereby vi- 
olated (B/C : model too low, 10 Be/ 9 Be: model too high[3 , local 
electron spectrum : model too high) this shows the difficulty of 
constructing such a model consistently. There are of course other 
possibilities to approach this issue, but we restrict ourselves to 
this example here. 

Another alternative to obtain the electron spectrum by prop- 
agation is to invoke an upturn in D(E) at low energies, in- 
stead of the constant used normally. Fig [12] shows a model with 
D(E) oc E~° 5 for E < 4 GeV, but otherwise the same as the 
pure diffusion model with low-energy injection index 2 and halo 
height 4 kpc. It reproduces the synchrotron data, and does not 
require the high electron spectrum of the model with small halo 
height. However it will under-predict the B/C data at low ener- 
gies - but this depends strongly on solar modulation and hence 
is not so critical. It is a possibly a more plausible scenario than 
the previous one since it violates less constraints. The required 
D(E) is in fact sim ilar to that in the wave-damping model of 
Ptus kin et alJ (|2006j) and hence has also a plausible physical ba- 
sis. 

Other variations on the propagation could affect the electron 
spectrum, for example spatial variations in the diffusion coeffi- 
cient, anisotropic diffusion and convection, which have not been 
considered here. They will not however affect our general con- 
clusion of the need for a significant break in the electron injec- 
tion spectrum. This is because any model must be constrained by 
the CR nuclei secondary /primary data, and the resulting mod- 
ified propagation parameters will finally lead to a path-length 
distribution similar to our basic model, and hence similar elec- 
tron energy losses. 

Although we might expect a dependence of D(E) on the B- 
field, in these models the variation of total B-field is very small 
(see section 2.3), so that including such a dependence would 



B/C oc zi,/D(E), 10 Be/ 9 Be oc \[D(E)/Zh. 
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have no sigificant effect. However in future such a dependence 
on B and its topology could be included to make the models 
more physically realistic. 

High-energy electrons are expected to have spatial variations 
due to their rapid energy losses combined with the stochastic 
nature of the sources in spac e and time. A study of this effect 
(Strong & Moskalenko 2001) shows that significant variations 
start above 10 GeV, and become large only above 100 GeV, and 
hence will not affect the low-energy synchrotron from lower en- 
ergies which is the main focus of this work. In any case the syn- 
chrotron spectra we use here are integrations over large sky areas 
and long lines-of-sight, which will smooth out any variations and 
hence not affect our conclusions. 



7. Discussion 

Independent of the propagation model, the primary electron 
spectrum must turn over below a few GeV, with an ambient in- 
dex around 2. In the pure diffusion model this implies an injec- 
tion index 1.3-1.6. It cannot cutoff completely since secondary 
leptons only contribute about one third of the low-frequency syn- 
chrotron, and electrons below 1 GeV are anyway directly ob- 
served by spacecraft in the heliosphere. 

While the synchrotron spectra may be subject to zero-level 
and scale errors, and the effects of absorption and free-free emis- 
sion may affect the comparison with models, the determination 
of the spectral index has been performed in a robust way by 
many authors and this gives a rather tight constraint on the am- 
bient electron spectrum. Th e most recent determ ination of the 
spectral index 45-408 MHz (Guzman et al.ll201ll) gives 2.5-2.6, 
implying an ambient electron inde J^l of 2.0-2.2 for electrons 
below a few GeV. This completely excludes a continuation of 
the ambient electron index 3.0-3.2 measured by Fermi-LAT >7 
GeV) to energies below a few GeV. Since in this range the dif- 
fusion coefficient is constant in the pure diffusion model, only 
energy losses steepen the spectrum slightly since escape domi- 
nates. This therefore implies an electron injection index of < 2, 
consistent with the detailed analysis presented here. While the 
latter is a very simplified argument it supports the conclusions 
of the detailed calculations in a robust model-independent way. 
At the same time the synchrotron index > 1 GHz has been found 
by many authors to be near 3, fully consistent with the measured 
ambient electron spectrum above a few GeV (ambient index 3.0- 
3.2 giving synchrotron index 3.0-3.1), and an injection index 2.2 
steepened by 0.5 from synchrotron/IC losses and by 0.5 due to 
D(E). 

To be more precise on this point, we can use an analytical 
approximation to the propagated electron spectrum for a plane 
parallel source distribution with diff usion (no rea c celer ation) 
and energy losses g iven e.g. by Bul anov & Dogell (1 19741) and 
iDelahave et al.l (l2010h . The spectrum steepens by (6 + a - l)/2 
where dE/dt oc E a and D(E)oc E s . For synchrotron losses a — 2, 
so at high energies where 6 = 0.5 the steepening is 0.75, in ac- 
cordance with our injection spectrum 2.2 and ambient spectrum 
3, and at low energies 6-0 the steepening is 0.50, consistent 
with our injection index 1.3-1.6 and ambient spectrum 2. In re- 
ality the relations are more complex due to spatial dependence 
etc., but the analytical form reproduces the general behaviour for 
synchrotron/IC-type losses, and shows it does not depend criti- 
cally on the detailed GALPROP modelling. 



16 The synchrotron spectral index for an ambient electron spectrum 
with power-law index p is approximately /?=2+(p-l)/2. 



Again we note that the determination of the ambient elec- 
tron spectrum is independent of the way in which this spectrum 
was produced via injection and propagation. We use a particu- 
lar type of parameterized model to generate physically plausible 
spectra, but other models would lead to the same ambient spec- 
trum since the final criterion is consistency with synchrotron and 
direct measurements. 

The low-energy interstellar spectrum has consequences for 
solar modulation: it must have a smaller effect than predicted 
by the force-field approximation for cD=600-700 MV generally 
used for such solar-mininum data (see Section 5); using our 
synchrotron-based spectrum, <$> <200 MV gives a better fit to 
data <2 GeV. The new PAMELA data reinforces this conclu- 
sion. Since the force-field approximation is anyway known to be 
unreliable, and the values used are often based on an assumed 
interstellar spectrum, we do not pursue this further, but simply 
propose that the interstellar spectrum determined in this paper 
be used in future physical modelling of modulation. From our 
analysis we just suggest that the low-energy falloff in the di- 
rectly measured electrons, normally attributed mainly to mod- 
ulation, may instead reflect more the interstellar spectrum. We 
note that although synchrotron probes CR leptons but not nuclei, 
the improved understanding of modulation would be relevant to 
all species. 

It is intentionally beyond the scope of this paper to specu- 
late on the origin of the electron injection spectrum; we sim- 
ply present it as an observational result posing a challenge 
for cosmic-ray source models. However some obvious remarks 
are in order. The low-energy injection spectrum is less than 
found in SNR; the distribution of SNR radio spectral indices 
IDelahave et al.l (1201 Oh is very broad - index 2.2 to 2.8, with mean 
2.48 giving an electron index 1.4 to 2.6, with mean 1.96. They 
quote an SNR electron index 2.0+0.3. However this is for elec- 
trons inside the SNRs: the escaping flux may be harder if higher- 
energy electrons escape more easily as might be expected. At 
high energies our derived injection index 2.2 suggests a steep- 
ening relative to the spectrum inside the remnants, but at lower 
energies this may not hold. A recent study of escape of electrons 
from SNR is given by Ohira et al. 201 1 darXiv:1106.1810l) . and 
lOhira et all d2010h , which includes among many effects an es- 
cape time decreasing with energy, and hence a hardening of the 
escape spectrum. Our results will constrain such models, which 
predict a complex escape spectrum. 

Other sources of electrons may also contribute, for ex- 
ample pulsar wind n ebula (PWN) have radio indices 2-2.3, 
(Reynolds et al. 201 1) giving electron index 1-1.6, similar to our 
low-energy injection index. While PWN clearly cannot produce 
the high-energy spectrum, it is possible that they play a role 
at low energies. Pulsars produ ce electrons with i ndex 1.5 - 1.9 
(with a cutoff at several GeV) (Gra sso et al ] |20IT1) again similar 
to our low-energy index. The attempt to construct a model with 
various source types to reproduce the data is beyond the scope of 
the present work, and in fact it is not easy to imagine how hard 
spectrum low-energy sources plus steeper spectrum high-energy 
sources could be combined to produce the observed composite 
spectrum. 

Finally we compare the electron injection spectrum with that 
of nuclei, in particular protons. In the GALPROP models used 
here, the nuclei injection spectrum has index, above/below a 
break energy of 9 GeV/nucleon, of 1.8/2.25 for the plain dif- 
fusion model, and 1.98/2.42 for the reacceleration model. These 
values have been chosen to agree with CR data for the given 
propagation parameters derived from B/C etc. Thus the low- 
energy nuclei index is slighty larger than what we have deduced 
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for electrons (1.3-1.6 for the plain diffusion model, 1.6 for the 
reacceleration model). However the nuclei are strongly affected 
by solar modulation at low energies, and unlike the electrons 
there is no equivalent of the synchrotron tracer for the interstellar 
spectrum (pion-decay gamma rays may be constraining but are 
sensitive mainly to protons above a few GeV). So a nuclei injec- 
tion spectrum equal to that of electrons would be acceptable, if 
indeed that should be predicted by a theory of CR acceleration, 
but is not required by the data. 

8. Conclusions 

Our main conclusion is that the interstellar CR electron spectrum 
must turn over rather sharply below a few GeV. This result is in- 
dependent of how the spectrum is formed by injection and prop- 
agation. The low-energy falloff in the directly measured elec- 
trons, normally attributed just to modulation, may instead reflect 
mainly the interstellar spectrum. The (model-dependent) injec- 
tion index implied for the primary electrons is 1.3-1.6 below a 
few GeV, and 2.1-2.3 at higher energies. The standard reaccel- 
eration model is not consistent with the observed synchrotron 
spectrum, since the total from primary and secondary leptons 
exceeds the measured synchrotron at low frequencies. While not 
excluding reacceleration models, it does pose a challenge to be 
addressed. 

We show that it is still possible to obtain the ambient electron 
spectrum by propagation even for a more conventional injection 
index of 2, but at the expense of violating other constraints. A 
low-energy upturn in the diffusion coefficient is the most promis- 
ing model of this kind. 

Therefore combining synchrotron data with direct measure- 
ments of CR provide unique and essential constraints on the in- 
terstellar electron spectrum. These results have impli cations for 
inters tellar gamma rays especially at low energies dPorter et al.l 
2008), and also for high energies (see Strong! (1201 ll) for a recent 
review). Exploiting the complementary information on cosmic 
rays and synchrotron gives new constraints and has implications 
for gamma rays. This connection is especially relevant now in 
view of the ongoing PLANCK and Fermi missions, and in fu- 
ture new radio astronomy instruments like LOFAR and direct 
measurements by AMS-02. 

There are of course limits to what conclusions based on a 
model like GALPROP can achieve since the full complexity of 
the Galaxy and physical processes can never be reproduced. A 
key is the combination of constraints from many different data 
types ('multi-messenger'), but still it is hard to break the de- 
generacy between the source injection spectrum and propagation 
even with secondary/primary ratios etc. In future as CR sources 
are better understood via radio, X-ray and gamma-ray observa- 
tions we should be able to make more basic progress on this 
topic. 
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Appendix A: Radio spectral index data. 

Starting with lower frequencies, Roge r et all (1 19991) used the 22 
MHz and 408 MHz surveys to deriv e average high latitude (3 = 
2.47, with a variation of about 0.05. lRogers & Bowman! 12008) 
measured over a continuous band 100-200 MHz, p — 2.5 ± 0.1 
; combining with other surveys 150-408 MHz p = 2.52 ± 0.04 
at high latitudes. This paper contains a useful summary table of 
P back to 1962. iGuzman etafj (1201 lh use the recent 45 MHz 
all-sky map together with the 408 MHz Haslam map to derive 
P = 2.5 - 2.6 over most of the sky for these f r equenci es. 

At higher frequencies, iReich & Reichl (Il988bl) . separat- 
ing thermal and nonthermal components for 408- 1420 MHz 
found a non-thermal index p = 2.85 - 3.1. IReich et all 
(12004 give spectral indices for 45-408-1420-22800 MHz, and 
give a table of zero-level corrections. Spectral index maps. 

(408 - \A2QMHz) = 2.6 - 2.7 away from plane and loops. 
Giardino et ail d2002l) find 408-1420 M Hz P = 2.78 ± 0.17 
and 408-2326 MHz p = 2.75 ± 0.12. iPlatania et all (12003b 
from a full-sky analysis of 408, 1420 and 2326 MHz surve 
find p 



G 



2.695 with a dispersion of 0.12. Pla tania et al. (199 
using radiometers at a high-altitude sit e at 1400-7500 MHz: 
P = 2.81 ± 0.16. IZannonietal.l (120081) report TRIS absolute 
measurements at 0.6, 0.82, 2 . 5 GHz , drift scans at 6 — 42° (th is 
is paper I, II ^Gervasi et all (l2008bl) . Il NTartari et ail d2008l) 1. 
iGervasi et al.l ( 2008al) derive the extragalactic sou rce contribu- 
tion to the background from 151 to 8440 MHz. iTartari et aD 
(2008) using TRIS III find a synchrotron halo with y8=2.9-3.1 
600-820 MHz. This paper discusses zero errors in other surveys 
(their table 9: 150, 408, 820, 1420 MHz surveys). They find 
that p increases from 2.2 to 2.8 from 150 MHz to 1420 MHz. 
They give a spectrum in two directions (9h,42°, 10h,42°) cor- 
rected for zero level (their fig 7) which can serve as a standard: 
it shows large steepening from 150 to 1420 MHz). They give 
the variation of ,6(600-820 MHz) along 6 = 42" : it is mainly in 
the range 2.8 to 3.2(their fig 5). This paper contains an exten- 
sive discussion of spectral indices past and present, and its rela- 
tion to electron spectrum. Tgal=T-Toff=T-(Tex+Tcmb)+Tzero. 
From their Table 9, the correction to the 408 MHz Haslam survey 
is Tzero=+3.9; EX=2.65 CMB=2.82 gives Toff=1.57 compared 
to lReich & Reichl d!988al) (their table 7): Tzero=+2.1,Toff=3.7± 
0.85. 



iKogut et al.l d2011h (ARCADE2) at 3, 8 , 10 GHz : j3 sync = 
2.55 ± 0.03 using a 408 MHz template, but they use cosec(b) 
and CII correlations. Poles/coldest regions p sync - 2.57 + 0.03 
These values are significantly lower than those normally found 
(seeabove:2.7-3.1).[3 

Turning to the highest frequencies, Hinsha w et al.l (120071) us- 
ing WMAP 3- year data found p = 3.15 - 3.5 for 23 - 61 GHz. 
iDunklev et al.l d2009l) used WMAP 5-year polarized maps to 
derive synchrotron p = 3.02 + 0.04 (at WMAP frequencies). 
They show s kymaps of and find no latitude dependence unlike 
IKogut et all d2007l) (WMAP 3-year data) who found an increase 
from 3 .05 to 3.25 from the Galactic plane to the poles. Gold et al. 
( 2009) used WMAP 5-year data to derive synchrotron p = 3.15+ 
0.10 (roughly from their fig 16), including polarized-only anal- 
ysis. T he latitude profi l e show s a lower index for \b\ < 10°: 2.8, 
(unlike IDunklev et al.l (|2009) wh o found no latitud e variation 
on the same 5-year data, but like Kogu t et al.l d2007l) who found 
simila r variation slightly different valu es) on the 3-year WMAP 
data). iMiville-Deschenes et al.l (120081) give maps of 23 GHz in- 
tensity and p 408 MHz-23 GHz using WMAP polarized emis- 
sion maps, especially useful for comparing with synchrotron 
models since it selects out the (polarized) synchr otron emission 
in a model-independent way. iTestori et al.l (120081) give spectral 
index distributions using polarized intensity for the 1 .435 GHz 
survey and 22 GHz WMAP. It peaks at/? = 2.7 - 3.0, but ranges 
from 1.8 to 3.6 due to depolarizatio n in the plane and being near 
the noise level. IReich et all d2004l) discuss zero-level errors in 
WMAP, and conclude, using T-T plots, that P( 1420 MHz - 22.8 
GHz) decreases from 3.1 to 2.8 after correction, nearer to the 
408-1420 MHz value. 



17 In a related ARCADE2 study, iFixsen et al] d201ll) claim a 3-90 
GHz extragalactic background about a factor 6 higher than expected 
from radio sources at 1 GHz. Excess claimed to have /?=2.60 from 22 
MHz to 10 GHz. To model the Galactic emission they use a cosec(b) 
analysis, which would not be sensitive to a large halo ( see their dis- 
cussio n, correlation with CII). This was followed up in Seiffe rTet al.1 
d20Hl) . who attribute the excess to underestimated Galactic emis- 
sion or unaccounted radio sources, or some combination of both. 
IVernstrom et all 1 201 1) use radio source counts to show that the claimed 
extragalactic background is hard to explain as the sum of sources to 
the currently observed flux limits, and that a extra po pulation at lower 
fluxes would be required. See also Singal et al. (2010), who also relate 
the radio to the extragalactic X- and gamma-ray backgrounds. 
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A. W. Strong et aL: Cosmic-ray electron spectrum from synchrotron 



galdef ID 54_z04LMPD_gO_1 .6_withoutsecS 



s galdef ID 54_zO4LllPD_g0_1.0_withsecS 




10 10' 10 10* 10 

frequency, MHz 



galdef ID 54_z04LMPD_gu_1 .6_onlysecS 



10 10' 10 10* 10 

frequency, MHz 



galdef ID 54 z04LMPDjO 1 .6 withsecS 



10° 




10° 




10' 10* 10° 
frequency, MHz 

Fig. 3. Synchrotron spectral index for pure diffusion model with 
primary low-energy electron injection index 1.6. Synchrotron 
from primary electrons (upper), secondary leptons (middle) 
and total (lower). The experimental ranges are based on values 
from the literature as reviewed in the text, and are only intended 
to be indicative of the general tren d since the measur ements 
cover different sky areas. Data: red: Tartari et al. (2008); blue: 
iRogers & Bowman! d2008l): cy an: iRoger et al.1dl999l): black 



dashed: Giardino et al. (2002); black dotted: IPlatania et al 
dl998 | ) : green: IPlatania et al l (120031): cyan: iKogutetaf. 
( 2007| ; orange full jMiville-Deschenes et all d2008l): orang e 



_ _ig e 

dashed iGold et all (l2009h; or ange dotted: Dunklev et al. (2009); 
black full :Kog ut et al.1 d201 ll) . 




Iff* 10 5 10 6 10 7 
energy, MeV 



Fig. 4. Electron spectra for pure diffusion model, low-energy 
electron injection index 1.0,1.3,1.6, 1.8,2.0,2.5. Modulation 
0=0,200,400,600,800 MV. Data as in ngureQ] 
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galdef ID 54_z04LMPD_gO_1 .0_wllhsecS 
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galdef ID 54_z04LMPD_gO_1 .3_wrthsecS 



i6 galdef ID 54_z04LMPD_gO_1 .8_withsecS 
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galdef ID 54_z04LMPD_gO_2.0_withsecS 




galdef ID 54 Z04LMPD jO 1 .6 wilhsecS 




Fig. 5. Synchrotron spectra for pure diffusion model with low- Fig. 6. Synchrotron spectral index for pure diffusion model with 
energy electron injection index (left to right, top to bottom) 1 .0, low-energy electron injection index (left to right, top to bot- 
1.3, 1.6, 1.8, 2.0, 2.5. Including secondary leptons. Data as in torn) 1.0, 1.3, 1.6, 1.8, 2.0, 2.5. Including secondary leptons. 
figure[2] Experimental ranges are based on the references reviewed in 

Section 4.1, and are intended to be representative not exhaus- 
tive. Data as in figure [3] 
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10 -21 y — I — I — I — I — I — r Vnm l 

1 10 1 2 1 3 1 4 1 5 1 6 
frequency, MHz 

Fig. 7. Synchrotron spectra for pure diffusion model with sharp 
cutoff in primary electrons below 4 GeV; primary electrons only 
(upper), primary and secondary leptons (bottom). The contribu- 
tion from secondary leptons is shown in Fig 1. Data as in fig- 
ure [2] 



Fig. 8. Lepton spectra for diffusive reacceration model with 
primary low-energy electron injection index 1.6. Primary 
electrons (upper), secondary positrons (lower). Modulation 
0=0,200,400,600,800 MV. Data as in figure^ 
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galdef ID 54_z04LM_withsecS 




1 10 1 2 10 3 1 0" 10 s 10 s 

frequency, MHz 



galdef ID 54_z04LM_withsecS 




1 10 10 ! 10 3 10" 10 s 10 s 
frequency, MHz 



Fig. 9. Synchrotron spectra for diffusive reacceration 
model with primary low-energy electron injection index 
1.6. Synchrotron from primary electrons (upper), secondary 
leptons (middle) and total (lower). Data as in figure[2] 



Fig. 10. Synchrotron spectral index for diffusive reacceration 
model with primary low-energy electron injection index 1.6. 
Synchrotron from primary electrons (upper), secondary leptons 
(middle) and total (lower). Data as in figure[3] 
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A. W. Strong et al. 



Cosmic-ray electron spectrum from synchrotron 




galdef ID 54_z01 LMPD_gO_2.0_withsecS 




1<T 10 3 
frequency, MHz 
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galdef ID 54_z04LMPDmod_gO_2.0_withsecS 
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Fig. 11. Top: Electron spectrum for pure diffusion model with Fig. 12. Top: Electron spectrum for pure diffusion model with 

halo height 1 kpc, low-energy electron injection index 2.0. low-energy D{E) cc E~ 0i , halo height 4 kpc, low-energy elec- 

Modulation (5=0,200,400,600,800 MV. Data as in figure [TJ tron injection index 2.0. Modulation 0=0,200,400,600,800 MV. 

Centre and bottom : Corresponding synchrotron spectrum and Data as in figure Q] Centre and bottom : Corresponding syn- 

spectral index, plotted as in figure [5] with which this figure chrotron spectrum and spectral index, plotted as in figure|5] with 

should be compared. which this figure should be compared. 



